Finite-bias conductance anomalies at a singlet-triplet crossing 



Chiara Stevanato, Martin Leijnse, Karsten Flensberg and Jens Paaske 
Center for Quantum Devices, The Niels Bohr Institute, 
University of Copenhagen, 2100 Copenhagen 0, Denmark 

Quantum dots and single- molecule transistors may exhibit level crossings induced by tuning ex- 
ternal parameters such as magnetic field or gate voltage. For Coulomb blockaded devices, this shows 
up as an inelastic cotunneling threshold in the differential conductance, which can be tuned to zero 
at the crossing. Here we show that, in addition, level crossings can give rise to a nearly vertical 
step- edge, ridge or even a Fano-like ridge-valley feature in the differential conductance inside the 
relevant Coulomb diamond. We study a gate-tunable quasidegeneracy between singlet and triplet 
ground states, and demonstrate how these different shapes may result from a competition between 
nonequilibrium occupations and weak (spin-orbit) mixing of the states. Our results are shown to 
be in qualitative agreement with recent transport measurements on a Mn complex [E. A. Osorio et 
al., Nano Lett 10, 105 (2010)]. The effect remains entirely general and should be observable in a 
wide range of Coulomb blockaded devices. 

PACS numbers: 85.65.+h, 33.80.Be, 73.23.Hk, 75.50.Xx, 



I. INTRODUCTION. 

Low-temperature transport measurements have re- 
vealed level crossings in Coulomb blockaded quantum 
dots and single-molecule transistors, induced either by 
tuning gate voltage^' or external magnetic field™!. 
These show up in a bias spectroscopy as inelastic cotun- 
neling thresholds crossing zero, when tuning the control 
parameter towards the level crossing. If the two differ- 
ent ground states have the same elastic cotunneling am- 
plitudes, the linear conductance shows hardly any trace 
of the crossing, except perhaps for a peak right at the 
crossing if the degeneracy gives rise to Kondo-effectP^. 
In the case of a transition from a singlet to a triplet 
ground state, however, the elastic cotunneling amplitude 
can easily be different on the two sides , and t he triplet 
state can even give rise to Kondo effectP^^. The dif- 
ference in elastic cotunneling leads to a sudden change 
in the linear conductance at the crossing, which may be 
smeared by a small mixing of the two competing ground 
states. In the case of a singlet-triplet crossing, a mixing 
due to spin-orbit coupling leads to an avoided crossing, 
which can be observed as a lower bound on the inelastic 
cotunneling thresholcP^. 

Two recent experiments on widely different single- 
molecule transistors have shown gate-induced singlet- 
triplet crossings. These experiments revealed an extra 
cotunneling line in the differential conductance pinned 
to the crossing at zero bias and extending up to bias 
voltages much larger than the other excited spin states, 
thus demarcating the boundary between two different re- 
gions 34 . In (the supplementary material to) Ref. [4j it 
was pointed out that such a feature could easily arise if 
the device has sufficiently asymmetric coupling to source 
and drain. In this case, even at finite bias voltage the 
molecule is largely in equilibrium with the strongest cou- 
pled electrode. Like at zero bias, a difference in the elastic 
cotunneling amplitude may therefore show up as a step 
edge extending into the finite bias region. This scenario 




FIG. 1: Sketch of differential cotunneling conductance as a 
function of bias voltage, V, and control parameter, x. Thick 
blue (solid) lines indicate inelastic cotunneling thresholds with 
an avoided crossing near zero bias, caused by a finite mixing, 
A, of the two involved states. The red (dashed) line indicates 
the line along which the levels cross, including the possibility 
of a slight bias dependence of the level splitting. The two 
different sides of this line may have a difference in dl/dV due 
to a difference in elastic cotunneling amplitudes for the two 
states. 



is sketched in Fig. [T] with differential conductance as a 
function of bias voltage V and control parameter x, which 
in these two experiments is gate voltage. In both exper- 
iments this line has a finite slope, implying that the bias 
voltage somehow influences the tuning of the crossing as 
part of the total effective gating of the relevant states. In 
this case, the step edge is accompanied by a peak deriving 
from differentiating a slanted step in current. In Fig. [I] 
the difference in conductance is reflected by two different 
shades of blue, and a conductance peak may be observed 
along the border marked by a red (dashed) line. This 



behavior can indeed be seen from Fig. S6a in Ref. [3j but 
the corresponding feature observed in Ref. |4] revealed a 
surprising dip-peak Fano-like structure, which remained 
unexplained. 

In this paper, we show that the dip-peak lineshape ob- 
served in Ref. 4 can arise from a competition between 
nonequilibrium effects and spin-orbit mixing. In general, 
a mixing, A, of the two states can lead to an avoided cross- 
ing as indicated in Fig. [I] However, in the experiment!^ 
the spin-orbit coupling is weak enough that no avoided 
crossing is observed in the inelastic cotunneling, i.e., it is 
smaller than the tunnel-induced broadening of the cotun- 
neling lines. Since the line is only weakly bias dependent 
(nearly vertical), its detailed line shape still provides oth- 
erwise inaccessible information about the spin-orbit cou- 
pling. We shall restrict our attention to the singlet-triplet 
crossing observed in Ref. HJ but the basic mechanism is of 
general validity and this phenomenon could be observed 
in a variety of cotunnel junctions. 

The paper is organized as follows. In Sect. [XT] we in- 
troduce the double-dot model which will be considered 
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throughout the paper. In Sect. [Ill] we present the results 
of a perturbative generalized master equation calculation 
of transport as a function of gate and bias voltage. This 
calculation reveals the singlet-triplet boundary in finite- 
bias cotunneling, and shows that it should connect to 
kinks in the Coulomb diamond edges arising from sequen- 
tial tunneling. In Sect. |IV[ we focus on the interior of the 
Coulomb diamond and recalculate the inelastic cotunnel- 
ing spectrum within an effective exchange cotunneling 
model. The effective cotunneling model is first estab- 
lished from a Schrieffer- Wolff transformation of the orig- 
inal double-dot model, 
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are 



and the features of Sect, 
reproduced. Finally, we include a singlet-triplet mixing 
in the form of a spin-orbit induced Dzyaloshinskii-Moriya 
term, and demonstrate how this causes a nonequilibrium 
population inversion which in turn gives rise to the peak- 
dip structure observed experimentally in Ref. HJ 



II. DOUBLE DOT MODEL 
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FIG. 2: (Color online) Sketch of double dot model. Electrons 
on different dots interact via their charge (Coulomb repulsion, 
U12) as well as via their spin (exchange, J12). Tunneling be- 
tween electrode r and dot i takes place with amplitude U r . 
The gate electrode is not shown. 



The double dot system is described by Hdd-> where 
is the energy of the single-particle orbital on dot z, which 
has number operator rii = ^ a ni a = J2a C L C ^ f° r s P m 
projection a =t>i- Ui and U12 are the onsite and inter- 
dot Coulomb charging energies, respectively. is the 
spin on dot i (where = if dot i is empty or doubly 
occupied). There is an interaction between the spins on 
the different dots, which can be dominated by either su- 
perexchange (J12 > 0) or Hund's rule coupling (J12 < 0), 
and arises as an effective low-energy description of a more 
complicated multi-orbital double dot model, where or- 
bitals on different dots are partially overlapping 4 . B is 
the magnetic field, which we assume to be applied in the 
^-direction and therefore coupled to the z-projection of 
the dot spins, Sf. 

Hr describes non-interacting electrons in reservoir 
r = L, i?, with energy e r k and number operator n r k a — 
c\ ka c r kcr- Tunneling between the dots and the electrodes 
is described by Ht- For simplicity we assume that the 
tunnel amplitudes, t r ^, are independent of k and that the 
density of states in the electrodes, p r , is energy indepen- 
dent. Then also the tunnel rates, T r i = 2tt p r \t r i\ 2 , which 
set the inverse time-scale for resonant single electron tun- 
neling (SET), are energy-independent. 



As a simple model which can undergo a transition be- 
tween singlet and triplet groundstates, we consider an 
interacting double quantum dot, where the two dots are 
connected in parallel to a left (L) and a right (R) elec- 
trode and capacitively coupled to a gate electrode, see 
Fig. [2] The Hamiltonian is given by H = Hdd + Hr + 
Ht, where 

H DD = ^2 e icr n icr + ^2 U i n it n H + ^12^1^2 

ia i 

+ Ji 2 Si-S 2 + S(5f + S|), (1) 
Hr = ^e r / c n r / ccr , (2) 

rka 

H T = ^ triCrkad^ + h.C. (3) 
rkia 



III. FULL GATE AND BIAS SPECTROSCOPY 

We calculate the nonequilibrium reduced density ma- 
trix of the double quantum dot, and the current flow- 
ing between the electrodes, /. using a generalized master 
equation (GME) approac h 1Q I 11 1 , based on the real-time 
transport theory! 2 -!. This approach treats all the local 
interactions (^,^12,^12) non-perturbatively. The tun- 
neling between the dots and electrodes is treated pertur- 
batively up to fourth order in Ht, which includes all 
coherent one- and two-electron tunnel processes, such 
as SET , elastic and inelastic cotunneling^, pairtunnel- 
ing } 14 | 15 l and tunnel-induced level broadening and shift. 
The Kondo effect is not included at this level of per- 
turbation theory, which is therefore applicable only at 
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FIG. 3: (Color online) Results of GME calculation presented 
as conductance maps, \dI/dV\ plotted on a logarithmic color 
scale as a function of V and V g (the absolute value is plot- 
ted since negative values, occuring in certain regions because 
of the tunnel-rate asymmetry, are not well represented on 
a logarithmic scale). The red arrows point out the singlet- 
triplet line. Here the tunnel amplitudes are tm = 30, t\L — 
0.7, tiR — 0.8, tiL — 0.5 (arbitrary units), the Coulomb charg- 
ing energies are U\\ — U22 = 2000T, U12 = 400T, and we 
choose e 2 = ei + 2T. B = in (a) and B = 50T in (b). 



temperatures much larger than the Kondo temperature. 
In fact, for the theory to be applicable also in the SET 
regime, the thermal energy scale must dominate over the 
tunnel broadening, T ^> T ri . 

The results of the GME calculations are shown in 
Fig. [3j represented as two-dimensional conductance 
maps, with the absolute value of the differential conduc- 
tance, \dI/dV\, plotted on a color scale as a function of 
bias voltage, V, and gate voltage, V g . Because of ca- 
pacitive effects, the energy of a state with N = J2i n i 
electrons on the double dot is ex —NV g , where we have 
set the proportionality constant (gate coupling) to one for 
simplicity. We assume equal capacitive couplings to left 
and right electrodes and apply the bias voltage symmet- 
rically with chemical potentials \±l,r = =tV/2, in which 
case the bias voltage does not result in direct capacitive 
energy- shifts. However, as was shown in Ref. [4i ; J 12 may 
depend on both V and V g (for simplicity we assume a lin- 
ear dependence). Therefore the N = 2 groundstate will 
change from singlet to triplet along a line in the (V,V g ) 
plane. 

Figure [3ja) shows the conductance within a voltage 
range where N = 2 in equilibrium. Already the rather 
simple model Hamiltonian ([!]) has a complicated many- 
body eigenspectrum and gives rise to a rich conductance 
map, but we focus here on the features originating from 
the transition between singlet and triplet ground states. 
The strong (yellow) conductance lines with a large gate 
dependence mark the onset of SET involving the N = 1 
(left) and N — 3 (right) charge states. The crossing 
of such lines at zero bias voltage indicates a charge- 
degeneracy point. Inside the central region (Coulomb 
diamond), the charge is fixed to N = 2, and the current 



is dominated by coherent two-electron cotunneling pro- 
cesses, where an electron is effectively transferred from 
one electrode to the other, involving only virtual occupa- 
tion of either the TV = 1 or N = 3 charge states. A co- 
tunneling process which does not transfer any net energy 
to or from the double dot is called elastic and gives rise 
to a constant background conductance. In an inelastic 
cotunneling process^, in contrast, energy is transferred 
to the double dot, which becomes possible when that 
energy can be supplied by the bias voltage. This gives 
rise to conductance features whenever the bias voltage 
equals an excitation energy of the double dot system. In 
first approximation the conductance feature is a step, but 
the joint effects of nonequilibrium population of excited 
states and Kondo correlations may give rise to a more 
complicated lineshape^. 

The inelastic cotunneling step corresponding to excita- 
tions from the singlet ground state to the triplet excited 
states emerges from the left diamond edges in Fig. [3^a). 
Because of the gate dependence of the singlet-triplet 
splitting (due to the gate dependence of J12), the inelastic 
cotunneling step acquires a finite slope and crosses zero 
bias close to the center of the Coulomb diamond. At this 
point the double dot ground state switches from singlet 
to triplet. To the right of this crossing a cotunneling step 
again emerges, corresponding to the opposite excitation 
(from the triplet, which is now the ground state, to the 
singlet). Also emerging from the zero-bias crossing is a 
weak conductance peak (for V > 0) or conductance dip 
(for V < 0), marked with red arrows, which signals the 
finite-bias position of the singlet-triplet transition (the 
slope is a result of the combined V and V g dependencies 
of J12). This "singlet-triplet line", which was experi- 
mentally observed in Refs. 3 4, is thus a signature of the 
ground state transition as a function of the applied volt- 
ages and appears since the elastic cotunneling conduc- 
tance of the triplet is somewhat larger than that of the 
singlet (the lineshape, and the conditions under which it 
is visible, are discussed in detail in Sect. IV). When the 
singlet-triplet line intersects the edge of the Coulomb di- 
amond, the slope of the diamond changes because the 
ground state is changed into one with a different voltage 
dependence. 

Figure [3jb) shows the same conductance map, but un- 
der the influence of a magnetic field. Those SET reso- 
nances which correspond to excitations with a finite total 
spin are now split by the Zeeman effect. In addition, the 
Zeeman splitting of the triplet state results in a splitting 
of the corresponding inelastic cotunneling steps. When 
the ground state is a singlet, excitations to the three 
triplet states lead to three distinct gate-dependent co- 
tunneling steps. To the right of the ground state crossing 
point, where the ground state is the T_i triplet, excita- 
tions to the singlet give rise to a gate-dependent cotun- 
neling step and excitations to the T triplet result in a 
gate-independent step (excitations from the T_i to the 
X+i triplet are suppressed by spin-selection rules for co- 
tunneling, AS Z = 0, ±1). In an experiment, the Zeeman 
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splitting of the inelastic cotunneling steps provides an im- 
portant check of the spin assignment of the excitations. 

In both Figs. |3ja) and (b), weak conductance peaks 
with the same gate dependence as the Coulomb diamonds 
are seen inside the Coulomb blockade region, close to the 
diamond edges. These peaks result from SET transitions 
starting from excited states^, which can be populated in 
a nonequilibrium situation due to inelastic cotunneling 
processes. 
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IV. LINESHAPE OF THE SINGLET-TRIPLET 
TRANSITION 



In Fig. [3j the singlet-triplet line appeared as a peak 
or a dip in the conductance, for respectively positive or 
negative bias voltage. This was also the result of the cal- 
culations presented in the supplementary information to 
Ref. [4j but as pointed out already there, the experimen- 
tal data from Ref. |4[ reproduced here in Fig. 4(a)| do not 
adhere to this behavior. Instead, the experiment shows 
a conductance peak, for both positive and negative bias 
voltage, which then evolves into a peak-dip, for V > 0, or 
dip-peak, for V < 0, roughly when the magnitude of the 
bias voltage becomes larger than all the relevant inelas- 
tic cotunneling thresholds. In this Sect., we will show 
that these Fano-like lineshapes along the singlet-triplet 
boundary can be reproduced by introducing a small mix- 
ing of the competing singlet and triplet states. In con- 
trast to the scenario in Fig. [3] and the one discussed in 
the supplementary information to Ref. |4l this effect relies 
entirely on nonequilibrium pumping of excited states and 
therefore shows up only in devices with not too different 
couplings to source and drain electrodes. 
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A. Effective exchange cotunneling model 
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We now focus our attention to the region deep in- 
side the Coulomb diamond with N = 2, where real 
charge fluctuations (SET) are suppressed, and second 
order processes (cotunneling), involving virtual fluctu- 
ations to the nearby N = 3 and N = 1 charge states, 
dominate. If Ui ^> U\2, there is always one electron in 
each dot when N = 2. In this limit, the system de- 
scribed by the Hamiltonian is effectively a two-spin 
system, and it is natural to work in the singlet-triplet ba- 
sis {|T + i), |Tq), |T!_i), \S)}. In this basis, the double-dot 
spin Hamiltonian reads: 



(4) 



FIG. 4: (a) The lower density plot shows d 3 I/dV 3 vs. bias 
voltage, V, and gate voltage, V g , for the Mn complex mea- 
sured in Ref. [U An external magnetic field of 10 Tesla clearly 
reveals a singlet groundstate with an excited triplet towards 
the left, and a triplet groundstate with an excited singlet to- 
wards the right. The zoom-in shows the same data now plot- 
ted as dl/dV. A steep faint line is seen to be pinned to the 
groundstate crossing. Cuts at various bias voltages, marked 
by dashed lines, reveal a distinct lineshape as a function of 
V g developing at large applied bias: dip-peak at V < and 
peak- dip at V > 0. The two nearly vertical dashed lines indi- 
cate the range in V g used for each of the cuts. The vertical 
range is the same in all the cuts, and has been centered at the 
mean value within each cut. (b) Density plot of differential 
conductance as calculated below in Sect. |IVB| with cuts at 
bias voltages V/T = ±269, ±76,0 (parameters are chosen as 
in Fig.|6|. The peak-dip (dip-peak) structure from panel (a) 
is clearly reproduced. 



5 



with eigenenergies 





_ 3 


(5) 




-- 4J12, 




= -B+±J 12 , 


(6) 






(7) 




=B+\j 12 , 


(8) 



where we have omitted the constant term U12 + e i + e 2- 
As discussed above, J 12 is assumed to depend linearly on 
V and V g , and taking the singlet-triplet crossing point to 
be the zero of our gate voltage axis, we thus assume an 
exchange coupling of the form: 



J12 = -OiVa - 13V, 



(9) 



in terms of two dimensionless numbers, a and (3. 

The kinetic energy of the electrode-electrons is still de- 
scribed by Hr and, applying a Shrieffer- Wolff transfor- 
mation 18 to effectively eliminate the charge fluctuations, 
we arrive at the following exchange cotunneling Hamil- 
tonian: 



cot 



r,r'=R,L 
i=l,2 

k,k' 



Wj, rr > \ t 

aa' 1 ^ <J(j/ J C rk<j Cr ' k 'cr' • 



Trrrr' 



(10) 



where r is the vector of Pauli matrices. Written in the 
singlet-triplet basis, the local spin-operators, 

Si = ((S+ + Sr)/2, -i(St - Sr)/2, S?) , (11) 
take the following form: 

S+ = ^=(\T +1 )(T \ + \T )(T. 1 \ (12) 
±\T +1 ){S\ T \S)(T_ 1 |), 

S~ = ^=(\T )(T +1 \ + \T. 1 )(T \ (13) 
±\S){T +1 \ T \T. 1 )(S\), 

SI = 1 -(\T +1 )(T +1 \-\T_ 1 )(T_ 1 \ (14) 
T|5)<r |T|2b>(5|). 

The spin exchange, and the potential-scattering ampli- 
tudes are related to the tunneling amplitudes appearing 
in ([3]) as: 



1 



1 



E2 — E% E2 — E\ 
1 1 



E2 — E3 E2 — Ei 



(15) 
(16) 



where the energy denominators keep track of the energy 
cost for the virtual charge fluctuations to respectively the 
N = 3 and the N = 1 states. Assuming that we are deep 
inside the Coulomb diamond, we have omitted the slight 
differences in these denominators for the different spin- 
states. The potential scattering term vanishes close to 
the center of the diamond, where E% — E2 = E\ — E2. 
Since this gate voltage dependence has little influence 
on the line shape of the singlet-triplet crossing studied 
below, one may, without loss of generality, assume that 
Ji.rr' — Wi^ rr ' . This will be done in some of the examples 
to follow. 

The cotunneling current can now be found as^ 



/ = e£(r^-r^)P(T?). 



(17) 



The non-equilibrium occupation numbers for the four 
singlet-triplet basis states, P(rf), are found by solving 
the rate equations 

= £ (T V , V P( V ) - T W P( V ')) , (18) 

77' 



with the constraint that 

i = E p M- 



and with rates 



(19) 



(20) 



determined from Fermi's golden rule as 

r^= 7 ^As^n B (A £ ^). (21) 

Here n# denotes the Bose function, and we have intro- 
duced the energy differences 



and the tunneling probabilities 



(22) 



<>, =v E E 

*=1,2 \j=x,y,z 



J, 



where the matrix elements T\ , 
off from Eqs. (fl2|)-(p|). 



(23) 

(7/ I Sf 1 77) can be read 



If the device is coupled much stronger to one electrode 
than the other, the nonequilibrium occupations of ex- 
cited states are strongly suppressed, and the finite-bias 
conductance becomes a simple succession of steps at bias 
voltages corresponding to energies of the excited states. 
Already from the elastic cotunneling amplitude, i^ L , 
one can see that the conductance is larger for a triplet 
than for a singlet ground state 



=47rp 2 Vn B (V)\W RL \ 2 , 



RL 

T_iT_ 



--4irp 2 Vn B (V) (\J RL \ 2 + \W RL \ 2 



(24) 
(25) 
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FIG. 5: (a) Differential conductance dl/dV vs. V/T at 
Vg/T = -300 (blue upper line) and at V g /T = 200 (red 
lower line). The curves are offset for clarity. The tunneling 
amplitudes are very asymmetric, tm = 30,£il = 0.7, tm — 
0.8, t2L — 0.5 (arbitrary units). We have furthermore chosen 
B/T = 50, a = 0.5, and f3 = 0.1. (b) Occupation probabili- 
ties of the different spin states at fixed bias V/T — 250 with 
parameters as in (a). The state |+) refers to the lowest lying 
state, i.e. \S) on the l.h.s. and \T-\) on the r.h.s.. |— ) refers 
to the next lowest, i.e. \T-±) on the l.h.s. and \S) on the r.h.s. 



FIG. 6: (a) Differential conductance dl/dV vs. V/T at 
Vg/T = -300 (blue upper line) and at V g /T = 200 (red lower 
line) . The curves are offset for clarity. The tunnel amplitudes 
are now chosen more symmetric than in Fig. [5] with tm = 0.3, 
and a finite spin-orbit coupling, D/T — 2\J~2 has been in- 
cluded. All other parameters are the same as in Fig. [5] (b) 
Corresponding occupation probabilities of the different spin 
states at fixed bias V/T — 250. Notice the occupation in- 
version between the two states |— ) and |+) at their would-be 
degeneracy point at V g /T = —150. 



B. Including singlet-triplet mixing 



which amounts to a factor of two in difference between 
conductance on the singlet, or the triplet side if Ji^r' — 
Wi :rr > . If the singlet-triplet boundary is completely ver- 
tical in the (V g , V)-plane (J12 depends only on V g ), this 
difference in conductance only shows up as a change in 
contrast in the conductance map as indicated in Fig. [I] 
If, on the other hand, the singlet-triplet crossing has a 
slight bias-dependence, this already explains why Fig. [3] 
shows a peak for positive, and a dip for negative bias 
voltage. This is seen clearly in Fig. [5|a), which shows 
dl/dV as a function of V for two different values of V g , 
calculated from formula ( fl7| ). The singlet-triplet cross- 
ing gives rise to a peak at positive bias (blue line) and 
a dip at negative bias (red line). Figure p|b), shows the 
corresponding nonequilibrium occupation numbers, and 
since we have chosen a much stronger coupling to the 
right electrode, nonequilibrium effects are weak and one 
observes only slight changes in a small region about the 
crossing. 



The cotunneling spectroscopy of the Mn complex in 
Ref. |4| reveals a transition between the lowest lying and 
the hi ghest lying triplet states (cf. r.h.s. of lower panel 
of Fig. |4(a)' ). This implies a violation of the normal spin- 
selection rules for cotunneling, from which we infer that 
the relevant spin-states are not really true eigenstates of 
S z , but instead a mixture of these inflicted by spin-orbit 
coupling in the molecule. We now add this information 
to our model by including a simple Dzyaloshinskii-Moriya 
ternPEZI H DM = (Si x S 2 ) • D, which reads 



a dm 



D 

2V2 



|T +1 ) + |T_ 1 ))(5| + |5)((T +1 | + (T_ 1 | 



(26) 



for D _L B. This has the important effect of mixing the 
competing singlet and triplet states and will in general 
give rise to an avoided crossing as illustrated in Fig. [l] 
with A ~ D. If, however, the coupling strength D is 
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smaller than the bias width of the cotunneling peaks, 
this avoided crossing need not be observed. Interestingly, 
even such a small spin-orbit mixing may still have an ob- 
servable effect on the line shape of the rather steep finite- 
bias singlet-triplet boundary. For simplicity we only in- 
clude the part of Hdm which mixes the singlet and lowest 
triplet (T_i) states, which is the only relevant term close 
to the singlet-triplet line. 

Including Hdm, we may determine the new eigen- 
states at a given bias voltage, V, which still enters via 
Ji2(V g , V). Using the same procedure as outlined above, 
the cotunneling current is readily recalculated. The re- 
sult obtained for D/T = 2\[2 is shown in Fig. [HJa), re- 
vealing a pronounced peak-dip, and dip-peak structure 
at respectively positive and negative bias, in qualita- 
tive correspondence with the experimental observation. 
Whereas the peak (dip) in Fig. [5| with no spin-orbit 
coupling included, could only be obtained with a sub- 
stantial difference in the couplings to respectively source 
and drain electrodes, the peak-dip (dip-peak) structure 
in Fig. [6] only arises when the couplings to source and 
drain are not too different. The observed peak-dip struc- 
ture is thus a nonequilibrium effect arising from the bias 
dependent occupation numbers of the various spin-states, 
shown in Fig.[6jb) to exhibit a population inversion of the 
two competing ground states near the crossing, around 
Vg/T ~ 50 in this plot. Apart from the requirement of 
nearly equal tunnel couplings, this population inversion 
and the concomitant peak-dip stru cture does not rely on 
any fine-tuning of parameters. Fig. 4(b) shows this result 
in a full density plot of the nonlinear conductance with 
a number of gate voltage traces taken at different bias 
voltages, all plotted in the same style as the experimen- 



tal data in Fig. 4(a) 



Reading off the slope of the nearly vertical singlet- 
triplet boundary and the slope of the inelastic singlet- 
triplet cotunneling threshold in zero magnetic field, one 
can now readily determine the constants a and j3 from 
formula Doing this, we find that a ~ 0.002 and 

j3 « 0.020. In Fig. |4(a)| the width in gate voltage of 
the peak-dip structure is roughly 0.05 V, and using a 
this now translates into a corresponding change in J12 of 
the order of 0.1 meV, which must be an indirect mea- 
sure of the effective singlet-triplet mixing, D, in this Mn- 



complex. 



V. CONCLUSIONS 

Inspired by the experimental results for the conduc- 
tance map of a Mn complex reported in Ref. |4j we have 
studied the signatures in the nonlinear conductance near 
a finite-bias singlet-triplet transition. For asymmetrically 
coupled devices the nonequilibrium effects are weak. A 
mere difference in elastic cotunneling amplitudes for re- 
spectively singlet and triplet states can give rise to cur- 
rent plateaus of different heights and hence a conduc- 
tance peak along the transition line in the (V, Vg)-plane. 
For more symmetrically coupled devices, this effect will 
be washed out at sufficiently large bias voltages. In this 
case, however, a new possibility arises, if the two compet- 
ing states are mixed by e.g. spin-orbit coupling. In this 
case, the nonequilibrium occupations exhibit a momen- 
tary population inversion near the transition line, which 
in turn gives rise to a peak-dip anomaly in the conduc- 
tance. This provides an explanation for the observations 
in Ref. 4, and it provides an indirect method of gauging 
the size of the singlet-triplet mixing. 

This type of cotunneling signature should be rather 
general, and the considerations about the singlet-triplet 
transition studied here could readily be exported to en- 
tirely different Coulomb blockaded systems with level 
crossings induced by gate- volt age. If the competing 
states are sensitive to gate voltage, it is only natural to 
expect also a slight sensitivity to bias voltage and thus a 
non- vertical transition line as depicted in Fig. [I] As was 
the case for the Mn complex discussed here, the steepness 
of this line can help bring out otherwise hidden spectro- 
scopic details. 
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